Limit on the Temporal Variation of the Fine-Structure Constant 

Using Atomic Dysprosium 
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Over a period of eight months, we have monitored transition frequencies between nearly degener- 
ate, opposite-parity levels in two isotopes of atomic dysprosium (Dy). These transition frequencies 
are highly sensitive to temporal variation of the fine-structure constant (a) due to relativistic cor- 
rections of large and opposite sign for the opposite-parity levels. In this unique system, in contrast 
to atomic-clock comparisons, the difference of the electronic energies of the opposite-parity levels 
can be monitored directly utilizing a radio-frequency (rf) electric-dipole transition between them. 
Our measurements show that the frequency variation of the 3.1-MHz transition in ^^''Dy and the 
235-MHz transition in ^®^Dy are 9.0±6.7 Hz/yr and -0.6±6.5 Hz/yr, respectively. These results 
provide a value for the rate of fractional variation of a of (—2.7 ± 2.6) x 10~^^ yr~^ (1 o") without 
any assumptions on the constancy of other fundamental constants, indicating absence of significant 
variation at the present level of sensitivity. 

PACS numbers: 06.20.Jr, 32.30.Bv 



Modern theories attempting to unify gravitation with 
the other fundamental interactions have renewed interest 
in experimental searches for temporal and spatial varia- 
tion of fundamental constants. These theories allow for, 
or even predict, such variationjwhich would violate Ein- 
stein's Equivalence Principle [l|. Recently, evidence for 
variation of a over cosmological time scales was discov- 
ered in absorption spectra of light from quasars [3, 
The result in Ref. [§| corresponds to a/a ~ (6.40±1.35) x 
10~^^/yr assuming a linear shift over 10^° years. How- 
ever, more recent measurements present sensitive results 
that are consistent with no variation 0, Q ■ On a geo- 
logical time scale of 10^ years, a test for variation of a 
comes from analyses of fission products of a natural re- 
actor in Oklo (Gabon). There is a discrepancy between 
earlier analyses Q , which are consistent with no vari- 
ation, and the analysis in Ref. Q , where a variation was 
reported. However, a recent full Monte Carlo simulation 
of the reactor reports a result consistent with no vari- 
ation at the level of 1.2 x 10~^''/yr. 

Observational measurements are difficult to interpret 
due to numerous assumptions and uncontrollable system- 
atic uncertainties. Laboratory searches (see, for exam- 
ple, Refs. 11, 12, l3|), which probe variations on the 
time scale of years, are easier to interpret since they are 
repeatable, and systematic uncertainties can be studied 
by changing experimental conditions. The best labora- 
tory limit of \a/a\ < 1.2 x 10~^^/yr was obtained from 
a comparison of a Hg+ optical clock to a Cs frequency 
standard [l3|- This result assumes that other constants, 
in particular me/mp, do not vary in time [l3 |. Inter- 
estingly, evidence for variation of this parameter on as- 
tronomical time scale was recently reported in Ref. [l5| . 



The best limit that is independent of assumptions re- 
garding other constants was obtained by combining this 
result with a comparison of Yb+ optical clock and a 
Cs frequency standard |13j] . This measurement gives 
\a/a\ < 2.2 x 10'^^ /yr. 

In this Letter, we report a limit on the temporal varia- 
tion of a that is independent of other constants. This re- 
sult is based on observing, over a period of eight months, 
rf-transitions between nearly degenerate, opposite-parity 
levels in atomic dysprosium, as suggested in Ref. |ld |. 

The energy of an atomic level can be expressed as 



E = Eo + qiayal-l), 



(1) 



where Eq is the present-day energy, ao is the present- 
day value of the fine-structure constant, and g is a level 
dependent coefficient that determines the sensitivity to 
variations in a [l7l |. A theoretical calculation has found 
that the ^-coefficients for the two nearly degenerate, 
opposite-parity levels in Dy (Fig. [2 levels A and B) are 
large and of opposite sign [t^ . For the even-parity level 
(level A), QA/hc ~ 6 x 10'^ cm^^, while for the odd- 
parity level (level B), qs/hc ~ —24 x lO'^ cm^^. Since 
explicit uncertainties are not quoted for these numbers, 
they are assumed to be exact in the determination of the 
result [13. 

Another parametrization of frequency and sensitiv- 
ity factor is often used for atomic-clock comparisons, in 
which an electronic transition frequency is expressed as 



RooAF{a), 



(2) 



where Roo is the Rydberg constant expressed in units of 
frequency, A is a dimensionless atomic structure factor 
independent of a, and F{a) is a factor that includes all 
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FIG. 1: Relevant levels and transitions in atomic dysprosium 
used to populate level B and detect the population of level 
A. Solid lines indicate laser excitation; dashed lines indicate 
spontaneous emission; wavy line indicates rf electric field. 



relativistic as well as many-body effects which depend 
on a [igj . The change in a transition frequency can be 
related to a fractional change in a as 



ly d 
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V at 
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a 



(3) 



The expression a d\nF/da is the sensitivity of a tran- 
sition to variation in a. Using this parametrization, the 
q-coefficients for atomic Dy can be recast (for a ~ ao) 
to give a d\nF/da = 2q/{hv), where hv is the energy 
of level A 01 B with respect to the ground state. This 
gives a sensitivity of 0.6 and —2.4 for level A and B, re- 
spectively. These values are comparable, for example, to 
the sensitivity of -3.2 for the transition used in the Hg+ 
optical clock [20| . 

In atomic-clock comparisons, the observable quantity 
is the ratio of the two frequencies being compared. A 
unique aspect of the Dy system is that the direct observ- 
able quantity is the difference vb — va, due to the fact 
that an electric-dipole transition can be induced between 
the two levels using a rf electric field. The time variation 
of the transition frequency between levels A and B is 



a 



-1.8 X lO^^Hz-. (4) 

a 



For instance, \a/a\ = 10~^^/yr implies |Ai^| ~ 2 Hz/yr. 

There are several advantages to using the nearly de- 
generate levels in Dy. Because most of the transition fre- 
quencies are ~1 GHz or smaller, direct frequency count- 
ing techniques can be used. This allows for the com- 
parison of two electronic transitions without the need 
for optical frequency combs or transfer cavities. A pre- 
liminary analysis of statistical and systematic uncertain- 
ties shows that the measurement of the transition fre- 
quency and the control of the systematics at a mHz level 



is feasible, which corresponds to an ultimate sensitivity 
of \ct/a\ ~ 10~^^/yr for two measurements separated by 
a year's time [2lj. A mHz resolution on a transition fre- 
quency of 1 GHz requires a relatively modest fractional 
stability of 10~^^ for the reference frequency standard. 
This also means that the results are insensitive to varia- 
tion of the Cs reference frequency due to changes in the 
values of fundamental constants since the experimental 
upper limit on thisj^lative variation rate is ~ 10^^^ yr^^ 
(see, for example, [ll]). In addition, because the isotope 
shifts and hyperfine splittings are on the order of the 
transition frequencies, it is possible to work with multiple 
rf transitions corresponding to the same electronic tran- 
sition, including, in particular, transitions with energy 
differences of opposite sign. Since the variation in a only 
depends upon electronic energy levels [l^], the variation 
of these rf transitions should have equal magnitude but 
opposite sign. This correlation can be used to detect and 
eliminate certain systematics. Currently, we monitor two 
such transitions: the 3.1-MHz {F = 10.5 F = 10.5) 
transition in ^^^Dy and the 235-MHz ( J = 10 ^ J = 10) 
transition in -'^^^Dy. 

The population and detection scheme is shown in 
Fig. [TJ The long-lived [tb > 200 /xs (H) level B is 
populated by three transitions. The first two transitions 
are driven with 833- and 669-nm laser light, while the 
third transition involves spontaneous decay. Atoms are 
transferred to level A with rf electric field referenced to 
a Cs frequency standard. Level A decays to the ground 
state in two steps. Fluorescent light at 564 nm from the 
second decay is used for detection. 

The atomic beam is described in Ref. [23|. It is pro- 
duced by an effusive oven operating at ~1500 K. In ad- 
dition to a multislit nozzle-array attached to the oven, 
two external collimators are used to coUimate both the 
atomic beam and the oven light (in order to minimize the 
background due to scattered light from the oven). The 
resultant atomic beam has a mean velocity of 500 m/s 
and a full-angle divergence of ^0.2 rad (1/e^-level). 

The laser light at 833 and 669 nm is produced by a 
Ti:Sapphire ring laser and a ring dye laser with DCM 
dye, respectively. Since narrow-band continuous-wave 
lasers are inefficient at exciting atomic beams with weak 
coUimation, an adiabatic passage technique is utilized to 
transfer the population to level / (a description of this 
technique, as well as references to earlier work, are given 
in Ref. [25|). Briefly, the beams are diverged with cylin- 
drical lenses to match the atomic-beam divergence. Due 
to the Doppler effect, the atoms experience a frequency 
chirp in laser detuning which transfers the population 
to the excited state with high efficiency. The laser-light 
and rf-interaction region is enclosed by a magnetic shield. 
The residual field in the interaction region is ~1 mG. 

The rf-generation and detection system is discussed in 
Ref. [2^. The frequency modulated rf field is generated 
by a synthesizer referenced to a commercial Cs frequency 
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standard which is compared to a second clock (Rb oscilla- 
tor locked to a digital cellular network signal) to monitor 
its stability. The short term fractional stability of the 
clocks is ~ 10~^^, and it approaches 10~^^ in 10 minutes 
of integration. The modulation frequency is 10 kHz with 
a modulation index of 1 . Since the rf-transition linewidth 
is ^^20 kHz, this modulation provides a fast sweep across 
the line shape, minimizing the effect of fluctuations such 
as those due to laser power drifts or density fluctuations 
of the atomic beam. The modulation is provided by the 
reference output of a lock-in amplifier which demodulates 
the signal from a photomultiplier tube used to detect the 
fluorescence. The line shape of the first- and second- 
harmonic outputs of the lock-in amplifier is, respectively, 
an odd function with zero crossing on resonance and an 
even function with a maximum on resonance. In order 
to reduce drifts further, the ratio of these two is used to 
measure the transition frequency, which is extracted by 
a two-step process described in detail in Ref . [26| . 

There are several imperfections that affect the stability 
of our transition- frequency measurements. The largest 
effect is due to the combination of the residual mag- 
netic field with the laser-light polarization imperfections. 
During normal operation, the polarizations of both laser 
beams are set to linear, which should lead to symmet- 
rically populated magnetic sublevels. However, residual 
light ellipticity leads to atomic orientation, which, in the 
presence of the residual magnetic field, leads to asymmet- 
ric broadening of the line and causes apparent shifts 21 1 . 
In systematic studies where we deliberately used circu- 
larly polarized light, shifts in the transition frequencies as 
large as ±160 Hz were observed. The linear polarizations 
of the laser beams are determined by calcite polarizers 
with 10~^ extinction ratio outside the vacuum chamber. 
However, the beams must travel through optical windows 
and lenses utilized for adiabatic passage population tech- 
nique (see above) before reaching the interaction region. 
A small amount of stress-induced birefringence on these 
lenses was discovered by deliberate misalignment of the 
laser beams. The amount of ellipticity depended upon 
where the beam sampled a lens. The residual systematic 
uncertainties due to this effect range from 2 to 5 Hz on 
different experimental runs with the smaller uncertainty 
corresponding to later runs in which the lens mounts were 
modified to relieve the stress. 

In addition to the residual magnetic field, the mag- 
netic shielding was inadvertently magnetized in the April 
2006 runs during studies utilizing external magnetic-field 
coils. Since the amount of magnetization and the resul- 
tant shifts in the transition frequencies were unknown at 
the time, systematic uncertainties have been assigned by 
driving the magnetization to saturation and noting the 
maximum shift induced for both field polarities. The un- 
certainty due to this effect is 7.5 Hz. The shields were 
demagnetized at the beginning of each subsequent run. 

Another important systematic uncertainty is due to 



the collisions of Dy atoms with the background gas in the 
vacuum chamber. CoUisional shift rates for the rf tran- 
sitions due to various gases were determined in Ref. [3] ■ 
The pressure in the vacuum chamber is ^ 5 x 10"'' Torr 
when the Dy oven is off. When the oven is turned on, 
the pressure rises to ^ 2 x 10^^ Torr, limited by H2 out- 
gassing from the oven. In Ref. [2^, the observed shift 
rates for H2 were found to be consistent with zero. How- 
ever, the pressure of other gases, such as N2, also in- 
creases at the level of~4xlO~^ Torr when the oven is 
turned on, and continually decreases and stabilizes after 
~ 6 hours of continuous oven operation. The partial pres- 
sures are monitored with a residual gas analyzer (RGA) , 
and the transition frequencies are corrected for the pres- 
ence of N2, O2, H2, and Ar using the shift rates measured 
for each gas. The variation in the pressure of other gases 
such as H2O and CO2 are ~ lO"'^ Torr or smaller, and 
no significant correlation between these pressure varia- 
tions and the transition frequency is observed. The total 
correction is < 2 Hz at any time during each run. The 
results also include small uncertainties (< 0.6 Hz) due to 
oven temperature variations, which may indicate a pos- 
sible effect of Dy intra-beam collisions [26| . 

Finally, there are systematic effects associated with rf 
electric-field inhomogeneities. The rf electrodes are fed 
by twisted-pair wires at one corner. Numerical simula- 
tions show that with this feed, at rf wavelengths compa- 
rable to the size of the electrodes (~ 10 cm), spatially 
varying amplitude and phase inhomogeneities become 
substantial. These inhomogeneities can lead to drifts in 
the transition frequency when they are combined with 
changes in the atomic beam velocity distribution. This 
effect was studied experimentally by deliberate detun- 
ing of the laser frequencies from the optical resonances. 
This leads to changes in the transverse velocity distri- 
bution of the excited atoms due to imperfect adiabatic 
passage. For the two transitions considered in this study, 
shifts of < 0.5 Hz for 1 MHz detuning in the laser fre- 
quencies were observed. To keep this effect at the 1 Hz 
level during the runs, the laser frequencies were periodi- 
cally retuned (~ every 2 minutes) to the center of their 
respective transitions. To check for slow variations in 
the beam velocity distribution due to other changes such 
as collimator clogging, this study was repeated several 
times over a two-year period, during which the oven was 
reloaded and collimators cleaned several times. No sig- 
nificant change in the size or dependence of the effect was 
observed. 

All other systematics considered in Ref. [213 were esti- 
mated to be smaller than 1 Hz for the specific conditions 
of the current experiment. Many of these systematics, 
such as AC Stark shifts and oven black-body radiation 
shifts, were studied experimentally by exaggerating cer- 
tain imperfections and were found to be negligible. 

We have measured the two rf transition frequencies 
in the course of seven runs over eight months. Fig- 
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FIG. 2: Measured transition frequencies for the 235-MHz and 
3.1-MHz transitions over an eight-month period. The data 
have been corrected for coUisional shifts. The soHd hnes are 
the least-squares linear fit to the data. The dashed lines are 
the least-squares fit to a constant function. 



ure[2]shows the results of these measurements corrected 
for coUisional shifts. Uncertainties are due both to 
statistical and aforementioned systematic uncertainties. 
Least-squares linear fit to the data points gives slopes 
of -0.6 ± 6.5 Hz/yr and 9.0 ± 6.7 Hz/yr for the 235- 
MHz and 3.1-MHz transitions, respectively. Accord- 
ing to Eq. (j4|), these results correspond to a /a = 
(-0.3 ± 3.6) X 10~^^ yr^i for the 235-MHz transition 
and (-5.0 ± 3.7) x lO'^^ yr'^ for the 3.1-MHz tran- 
sition. As mentioned earlier, because the energy dif- 
ferences for these two transitions are of opposite sign, 
these two data sets can be combined to construct sum 
and difference frequency plots. The sum frequency 
slope, which should be insensitive to a variation in a, 
is 8.1 ± 9.4 Hz/yr and consistent with zero. The differ- 
ence frequency, which is twice as sensitive to a variation 
in a as the individual frequencies, gives the final result of 
a/a = (-2.7 ± 2.6) x 10"^^ yr'^ consistent (1 a) with 
no variation at the present level of sensitivity. 

In conclusion, we have presented the first result of a 
direct measurement of the temporal variation of a with 
atomic dysprosium. This result is of comparable un- 
certainty to the present best laboratory result that is 
independent of other fundamental constants. However, 
in our case the interpretation does not require compari- 



son with different measurements to eliminate dependence 
on other constants, and provides an alternative to mea- 
surements that utilize state-of-the-art atomic optical fre- 
quency clocks. The present uncertainty is dominated 
by systematic effects primarily due to polarization im- 
perfections coupled to the residual magnetic field, col- 
lisional shifts, and rf electric-field inhomogeneities. A 
significant improvement is expected from a new appara- 
tus under construction which will provide better control 
over these, as well as other, systematic effects expected 
to be important to achieve better than 1-Hz sensitivity. 
Ultimately, mHz-level sensitivity may be achievable with 
this method 21 1. 
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